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ABSTRACT 

Intervening CIv absorbers are key tracers of metal-enriched gas in galaxy haloes over cosmic time. Previous studies suggest 
that the CIV cosmic mass density (Qc1y) decreases slowly over 1.5 S z <5 before declining rapidly at z = 5, but the cause of 
this downturn is poorly understood. We characterize the cy evolution over 4.3 S z S 6.3 using 260 absorbers found in 42 
XSHOOTER spectra of z ~ 6 quasars, of which 30 come from the ESO Large Program XQR-30. The large sample enables us to 
robustly constrain the rate and timing of the downturn. We find that cy decreases by a factor of 4.8 + 2.0 over the ~300 Myr 
interval between z ~ 4.7 and ~5.8. The slope of the column density (log N) distribution function does not change, suggesting 
that C IV absorption is suppressed approximately uniformly across 13.2 < log N/cm? < 15.0. Assuming that the carbon content 
of galaxy haloes evolves as the integral of the cosmic star formation rate density (with some delay due to stellar lifetimes and 
outflow travel times), we show that chemical evolution alone could plausibly explain the fast decline in Qc, over 4.3 S z < 6.3. 
However, the C IV/C II ratio decreases at the highest redshifts, so the accelerated decline in Qc y at z 2 5 may be more naturally 
explained by rapid changes in the gas ionization state driven by evolution of the UV background towards the end of hydrogen 
reionization. 
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observational claims of unexpectedly massive galaxies at z = 10 


E ERO revealed by JWST (e.g. Atek et al. 2023; Furtak et al. 2022; Harikane 


The formation of the first galaxies marked an important turning 
point in cosmic history. Massive stars released high energy photons 
which commenced the reionization of the Universe, and stellar 
nucleosynthesis led to the production of the first heavy elements, 
which were then released into the surrounding gas via supernova 
explosions. However, relatively little is known about the timing 
of the formation of the first galaxies and how they shaped the 
properties of their surrounding environments (see Bromm & Yoshida 
2011; Dayal & Ferrara 2018; Robertson 2022, for reviews). The 
shortcomings in our understanding have been emphasized by recent 
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et al. 2022a; Labbe et al. 2022; Naidu et al. 2022; Yan et al. 2022; 
Adams et al. 2023). 

Our understanding of the properties of the Universe near the end 
of the Epoch of Reionization has been propelled by the detection 
of growing numbers of z ~ 6 quasars. Constraints on the timing 
of reionization from measurements of Lya forest dark gaps (e.g. 
Zhu et al. 2021; Zhu et al. 2022) and transmission statistics (e.g. 
Fan et al. 2006; Becker et al. 2015b; Bosman et al. 2018; Eilers, 
Davies & Hennawi 2018; Yang et al. 2020; Bosman et al. 2022), and 
Lya damping wing absorption (e.g. Davies et al. 2018; Greig et al. 
2022) are consistent with leading independent probes of reionization 
(e.g. the luminosity function of Lya emitters; Konno et al. 2018, 
the fraction of Lyman Break galaxies showing Lyo emission; Mason 
et al. 2018, and the optical depth to reionization; Planck Collaboration 
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et al. 2020), and suggest that signatures of reionization persist below z 
~ 6. Furthermore, intervening metal absorption lines in z ~ 6 quasar 
spectra provide key insights into the chemical content and ionization 
state of gas around early galaxies. The optical depths of absorption 
lines remain constant as the quasar light travels through space, 
making it possible to investigate the properties of faint low-mass 
galaxies that fall below the detection limits of current emission-line 
surveys. The comoving mass densities of metal ions trace the overall 
chemical content of the Universe over cosmic time, whilst ionic ratios 
provide constraints on the ionization state of the absorbing gas (see 
e.g. Becker, Bolton & Lidz 2015a for a review). 

The C1v 441548, 1550 Å doublet is commonly used as a tracer 
of enriched gas across cosmic time because it is observable across a 
wide redshift range and is easily identifiable due to its doublet nature. 
CIV is produced by photons with an energy of at least 47.9 eV and 
primarily traces metals in the circumgalactic and intergalactic media 
(CGM and IGM), observed as absorption systems (e.g. Pettini et al. 
2003; Schaye et al. 2003; Péroux et al. 2004; Simcoe, Sargent & 
Rauch 2004; Songaila 2005, 2006; Simcoe 2006; Schaye, Carswell & 
Kim 2007; Danforth & Shull 2008; Becker, Rauch & Sargent 2009; 
Ryan-Weber et al. 2009; D'Odorico et al. 2010; Simcoe 2011; Tilton 
et al. 2012; Cooksey et al. 2013; D’Odorico et al. 2013; Shull, 
Danforth & Tilton 2014; Boksenberg & Sargent 2015; Burchett et al. 
2015; Danforth et al. 2016; Diaz et al. 2016; D’Odorico et al. 2016; 
Kim, Carswell & Ranquist 2016; Bosman et al. 2017; Codoreanu 
et al. 2018; Cooper et al. 2019; Hasan et al. 2020; Manuwal et al. 
2021; Hasan et al. 2022) as well as emission nebulae around quasars 
(e.g. Guo et al. 2020; Travascio et al. 2020). 

Leading theories of CGM and IGM enrichment suggest that the 
majority of the carbon traced by CIV absorbers was not formed 
in situ, but was produced in stars and subsequently ejected from 
galaxies by means of outflows (e.g. Aguirre et al. 2001; Theuns et al. 
2002; Oppenheimer & Davé 2006; Kobayashi, Springel & White 
2007; Oppenheimer, Davé & Finlator 2009; Cen & Chisari 2011; 
Finlator et al. 2013, 2020; Yamaguchi, Furlanetto & Trapp 2022). At 
z 25, star-formation driven outflows are expected to be the dominant 
source of metal-enrichment (e.g. Tescari et al. 2011; Pallottini et al. 
2014; Suresh et al. 2015; Sorini, Davé & Anglés-Alcázar 2020). 
Measurements of the statistics of CIV absorbers at z = 5 therefore 
provide valuable constraints on subgrid feedback models, which are 
typically tuned to reproduce the star formation rate and stellar mass 
content of galaxies but not the CGM properties (e.g. Oppenheimer & 
Davé 2006; Suresh et al. 2015; Keating et al. 2016; Rahmati et al. 
2016; García et al. 2017b; Finlator et al. 2020). 

Previous studies of C Iv absorption across cosmic time have found 
that the comoving mass density of CIV (Qc) declines smoothly at 
1.5 S z S 5 before dropping rapidly at z 2 5 (e.g. Songaila 1997, 
2001, 2005; Becker et al. 2009; Ryan-Weber et al. 2009; D'Odorico 
et al. 2010; Simcoe 2011; D'Odorico et al. 2013; Boksenberg & 
Sargent 2015; Díaz et al. 2016; Bosman et al. 2017; Codoreanu et al. 
2018; Meyer et al. 2019; D'Odorico et al. 2022). This could be an 
indication that the carbon content of galaxy haloes grows quickly 
over 6 = z 2 5 due to early enrichment by outflows. However, 
Qcıy also depends on the fraction of carbon existing as CIV, which 
is determined by the ionization state of the CGM/IGM. There is 
growing observational evidence for a rapid transition in the typical 
ionization environments of metal absorbers at z ~ 5.7. The observed 
decline in C Iv absorption is mirrored in SiIv 441393, 1402 A (e.g. 
D'Odorico et al. 2022), which is produced by photons of similar 
energy to CIV (733.5 eV). In contrast, the incidence of weak low- 
ionization absorbers (e.g. Mg Il 142796, 2803 A and O 1 41302À) 
remains constant or increases at z = 5.7 (e.g. Bosman et al. 2017; 
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Chen et al. 2017; Codoreanu et al. 2017; Becker et al. 2019). The 
ratio of CIV to C11 11334 A, which probes the distribution of carbon 
between ionization states, drops at z = 5.7 (e.g. Cooper et al. 2019). 
CTV absorption is often weak or undetected in low-ionization metal 
absorbers at z ~ 6 (e.g. D'Odorico et al. 2018; Simcoe et al. 
2020), unlike similar systems at z ^ 2-3 which show ubiquitous 
Civ absorption (e.g. Songaila 2006; Prochaska, Hennawi & Simcoe 
2013; Rubin et al. 2015; Cooper et al. 2019). These findings suggest 
that there is a transition in the average ionization environments of 
metal absorbers, from more neutral and low-ionization environments 
at z ~ 6 to more highly ionized environments at z S 5. This 
transition may be driven by an increase in the hardness and/or 
amplitude of the UV background following the end of reionization 
(e.g. Oppenheimer et al. 2009; Finlator et al. 2015; Becker et al 
. 2019). 

Large quasar samples are required to robustly examine the relative 
contributions of chemical enrichment and changes in ionization state 
to the rapid evolution in Qc over 5 S z S 6. The Ultimate 
XSHOOTER legacy survey of quasars at z ^ 5.8-6.6 (XQR-30; 
D'Odorico et al. in prep) has obtained 30 high signal-to-noise 
(S/N) quasar spectra at ^-30 km s^! resolution, nearly quadrupling 
the previous sample of 12 high quality z ~ 6 quasar spectra in 
the XSHOOTER archive. The combined sample of 42 spectra is 
referred to as the enlarged XQR-30 or E-XQR-30 sample. The 
first publications based on this sample have revealed that relativistic 
quasar-driven outflows are much more prevalent at z ~ 6 than at z 
~ 2-4 (Bischetti et al. 2022), that fluctuations related to re-ionization 
persist in the IGM until at least z — 5.3 (Zhu et al. 2021; Bosman 
et al. 2022), that bright z ^ 6 quasars typically live in overdense 
environments (Chen et al. 2022), and that quasar broad line regions 
atz ~ 6 are already enriched to ~2—4 times the solar abundance (Lai 
et al. 2022). 

Davies et al. (2023, hereafter Paper I) performed a systematic 
search for metal absorption lines in the E-XQR-30 spectra and 
published a catalogue of identified absorption systems. We use a 
sample of 260 intervening C Iv absorbers from this catalogue to make 
robust measurements of the number density, cosmic mass density, and 
column density distribution function of C Iv absorbers over 4.3 S z S 
6.3. Using these measurements, we examine when and how rapidly 
the drop in Qc, occurs, and explore whether this decline is driven 
primarily by chemical evolution and/or a change in the fraction of 
carbon found as C Iv. 

The paper is structured as follows. We describe the metal ab- 
sorber catalogue and the properties of the CIV absorber sample in 
Section 2.3. The CTV line statistics are presented and compared 
to predictions from the Technicolor Dawn simulation in Section 3. 
We examine the contributions of chemical evolution and changes in 
ionization state to the Qc,y evolution in Section 4, and summarize 
our results in Section 5. 

Throughout this work, we adopt the Planck Collaboration et al. 
(2020) ACDM cosmology with Ho = 67.7 km s^! Mpc^! and Qm 
031. 


2 SAMPLE AND DATA PROCESSING 


2.1 Quasar spectra 


The sample of C Iv absorbers used in this paper is drawn from the 
publicly available XQR-30 metal absorber catalogue,! which was 


! Available on GitHub: https://github.com/XQR-30/Metal-catalogue 
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assembled through a systematic search for absorption lines in the 
spectra of 42 z ~ 6 quasars (Paper I). 30 of the quasars were observed 
as part of XQR-30, an ESO Large Program (PI: V. D'Odorico) which 
obtained deep (continuum S/N ratio >10 per 10 km s^! spectral pixel 
at 1285 Å rest-frame), medium resolution (R ~ 10000, FWHM ~ 
30km s^), spectra of luminous (Jag < 19.8) quasars at 5.8 S z 
S 6.6 using XSHOOTER (Vernet et al. 2011). The observations 
were performed using slit widths of 0.9 and 0.6 arcsec for the VIS 
and NIR spectroscopic arms, respectively. The remaining 12 spectra 
were sourced from archival XSHOOTER observations of quasars in 
the same redshift and J-band magnitude range that were observed 
at similar spectral resolution and meet the S/N requirements for 
XQR-30. The full data set is described in detail in D'Odorico et al. 
(in prep). All quasars were selected based on their redshift and J- 
band magnitude, with no prior knowledge of intervening absorber 
properties. 

The XQR-30 and archival observations were reduced using the 
same pipeline for consistency. The data reduction process is outlined 
in Becker et al. (2019). Briefly, a composite dark frame was 
subtracted from each exposure, sky subtraction was performed on 
the un-rectified frame (Kelson 2003), and an initial ID spectrum was 
extracted using the optimal weighting method of Horne (1986). The 
Cerro Paranal Advanced Sky Model (Noll et al. 2012; Jones et al. 
2013) was used to compute telluric corrections which were then 
applied back to the 2D spectra. To maximize the bad pixel rejec- 
tion efficiency, all exposures were processed simultaneously when 
extracting the final 1D spectrum for each quasar and XSHOOTER 
spectroscopic arm (VIS and NIR). The final extracted spectra have a 
velocity sampling of 10 km s^!. 

The widths of the absorption line profiles revealed that when 
the average seeing is significantly smaller than the slit width, the 
true spectral resolution can be significantly higher than the nominal 
(slit-width-dependent) value. The true spectral resolution of each 
extracted spectrum was estimated as described in D'Odorico et al. 
(in prep). The average widths of spectral order spatial profiles were 
measured for individual high S/N 2D frames (a single value per 
exposure), and used to derive empirical relationships between the 
average spectral resolution and the seeing for each XSHOOTER 
arm and slit width. The empirical relationships were then used to 
convert the recorded seeing values into estimates of the Gaussian 
width of the line spread function (LSF) for each exposure. The final 
spectral resolution for each spectrum was obtained by computing 
the (S/N?-weighted average of the individual LSFs. The adopted 
spectral resolutions are listed in Paper I and range from R = 9500- 
13700 (median 11400, corresponding to FWHM = 26km s^!) for 
the VIS arm and R = 7600-11 000 (median 9750, corresponding to 
FWHM = 31 km s7!) for the NIR arm. 

Finally, the VIS and NIR spectra were combined to pro- 
duce a single spectrum for each quasar The NIR spectrum 
was scaled to match the flux level of the VIS spectrum 
over 990-1015 nm, and both spectra were cut at 1015nm be- 
fore being stitched together, manually correcting for small (Av 
~ 5-20kms ^!) mismatches between the relative wavelength 
calibrations of the arms as necessary We note that the ab- 
solute flux scaling of the combined spectrum does not im- 
pact the derived absorption line properties which are mea- 
sured from continuum-normalized spectra. The continuum fit- 
ting was performed using univariate spline interpolation within 
the PYTHON library ASTROCOOK? (Cupani et al. 2020), as de- 
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scribed in Paper I. The continuum fits are publicly available 
at the same location as the metal absorber catalogue (see 
footnote 1). 


2.2 Metal absorber catalogue 


The 42 quasar spectra were systematically searched for metal absorp- 
tion lines as described in Paper I. Briefly, an initial candidate list was 
generated by peforming an automated search for systems showing 
absorption in Mg 1, Fe I1 (4.4.2344, 2382, 2586, and/or 2600À), Cian 
conjunction with other low-ionization lines), C IV, Si Iv, and/or N V 
(44.1238, 1242 A) using ASTROCOOK. Only the wavelength regions 
redward of the quasar Lyo emission lines were searched because the 
saturation of the z ~ 6 Lyo forest makes it impossible to identify 
individual absorbers at shorter wavelengths. When performing this 
search, wavelength regions showing particularly strong skyline or 
telluric residuals in individual quasar spectra (most commonly the 
regions between the J, H, and K observing bands at 1.35-1.45 um and 
1.8-1.95 um) were masked to prevent the algorithm from returning a 
large number of spurious candidates. We note that approximately half 
of the spectra show broad absorption line (BAL) features associated 
with quasar-driven outflows (Bischetti et al. 2022 and submitted). 
The BAL regions were not masked, but any absorbers found in these 
regions were flagged in the final catalogue. 

The list of candidate systems was automatically filtered to remove 
spurious systems originating from strong skyline residuals or chance 
alignment of unassociated transitions in velocity space. For each of 
the remaining candidates, an automated search was performed to 
identify potential absorption in additional ions and/or transitions at 
the same redshift. Regions showing strong skyline/telluric residuals 
were not masked during this search because the redshift prior set 
by the initial system significantly increases the probability that any 
aligned absorption is real. The candidate systems were visually 
inspected by five experts to remove remaining spurious systems. 
After this initial line identification phase was completed, the spectra 
were searched for unidentified absorption features. In cases where 
unidentified absorption lines could be visually identified, the relevant 
transitions/systems were added to the catalogue with flags indicating 
that they were not recovered by the automatic line finder. 

The final metal absorber catalogue contains a total of 727 CIv 
absorption components spanning the redshift range 4.353 <z < 
6.487. The physical properties (z, Dopper b parameter, and column 
density, log N = log N/cm?) of all absorption components were 
measured by fitting Voigt profiles using the ASTROCOOK GUI. 


2.3 Sample selection 


In order to accurately measure the evolution in CIV absorber 
properties over cosmic time, we require robust estimates of the 
sample completeness (the fraction of true systems that are recovered) 
and the false-positive rate (the fraction of recovered systems that are 
spurious) as a function of log N. The calculations are described in 
detail in Paper I. A set of 840 mock spectra (20 per quasar) was 
generated with known absorber properties (spanning a wide range in 
redshift, log N and b parameter). The mock spectra were processed 
using the same steps applied to the observed spectra to generate mock 
absorber catalogues, which were then compared to the lists of inserted 
systems. The only exception is that no attempt was made to visually 
determine the origin of unidentified absorption features in the mock 
spectra. Therefore, the completeness and false-positive rate are only 
well defined for systems that were found in the automated line search 
and pass the visual inspection check. The completeness cannot be 
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Figure 1. Left: Illustration of the C Iv absorber sample. Horizontal lines show the redshift intervals over which the search for C Iv absorbers was conducted 
for each of the 42 quasar sightlines. Grey-dashed regions highlight proximity zones within 10000 km s^! of the quasar redshift. Longer gaps trace redshift 
intervals where the C Iv lines fall within BAL features or regions that were masked due to strong skyline or telluric contamination. The markers show all C Iv 
absorbers in the catalogue, where solid circles indicate primary absorbers (those that were automatically identified, pass the visual inspection check, and do not 
fall in masked wavelength regions or BAL regions) and open squares indicate secondary absorbers (all others). Proximate absorbers are shown in black. The 
intervening absorbers are split into two redshift bins (indicated by the marker colour), divided at the path-length-weighted mean redshift of our survey. Right: 
Completeness as a function of column density for C Iv absorbers in the two redshift intervals considered in Paper I. The dashed lines show the best-fit arctan 


functions published in that work. 


robustly quantified in BAL regions because the broad absorption 
makes it more difficult to recover underlying narrow absorption. The 
completeness in the masked wavelength regions is poorly defined 
because these regions were excluded from the initial line search. 

For these reasons, the metal absorber catalogue is divided into a 
primary sample which contains absorbers that were automatically 
identified, pass the visual inspection check, and do not fall in 
masked wavelength regions or BAL regions, and a secondary sample 
containing absorbers that fail one or more of these criteria. In this 
study, we only use the primary absorber sample to ensure that the 
completeness corrections are robust. 

In the left-hand panel of Fig. 1, the horizontal lines show the 
redshift ranges over which the search for CIV absorption was 
conducted for each of the 42 quasar sightlines. The grey-dashed 
regions highlight the proximity zones (within 10000 km s^! of the 
quasar redshifts), and longer gaps trace redshift intervals, where 
Civ falls within BAL features or regions that were masked due to 
strong skyline or telluric contamination. The markers show all of the 
CIV systems in the catalogue, where solid circles and open squares 
represent absorbers that are part of the primary and secondary sample, 
respectively. Of the 727 CIV components in the complete catalogue, 
559 are part of the primary sample. 

We additionally restrict our sample to intervening absorbers 
(which are conservatively defined to have a velocity offset of 
710000 km s^! from the quasar redshift) because absorbers in close 
proximity to quasars may not reflect the average ionization state and 
abundance patterns of the underlying absorber population (e.g. Berg 
et al. 2016; Perrotta et al. 2016). Proximate absorbers are shown 
in black, while the intervening absorbers are shown in blue and 


5We examined the impact of this choice on the derived C Iv absorber statistics 
by measuring the number densities and cosmic mass densities of strong C IV 
absorbers (13.8 < log N < 15.0, for which the completeness corrections 
are neglible) i) using only the primary sample, and ii) including manually 
identified systems that do not lie in masked/BAL regions. The two sets of 
measurements differ by — 15 per cent, and are consistent within the lo errors. 


red (where the colour indicates whether the absorber lies above or 
below the path-length-weighted mean redshift of the survey; see 
Section 2.6). Of the 559 primary C Iv absorbers, 507 are intervening. 

Paper I showed that the false-positive rate in the catalogue is 
negligible (<5 per cent) at all column densities, and the completeness 
as a function of log N is well described by an arctan function. The 
best-fit arctan parameters were provided in two redshift bins (4.32 
<z < 546 and 5.46 < z < 6.76) as shown in the right-hand 
panel of Fig. 1. The C v sample is 50 percent complete for log N 
> 13.08 (logN > 13.22) at 4.32 <z < 5.45 (5.45 < z < 6.76, 
see table 2 of Paper I). The survey is slightly less sensitive to CIV 
absorption in the higher redshift range because the typical noise 
level in the spectra increases towards longer wavelengths. In our 
analysis, we only include components with column densities above 
the 50 percent completeness limit of the higher redshift bin (log N 
> 13.2) because the number of detected systems drops at lower 
column densities, increasing the sampling error and reducing the 
reliability of the measurements. This cut reduces the sample from 
507 components to 274. 

Finally, we exclude components with log N > 15 to be consistent 
with previous works. This criterion results in the removal of a 
single component at z = 5.109 in VDESJ0224-4711. Our final 
sample consists of 273 CIV absorption components spanning the 
redshift range 4.353 <z < 6.067. Table 1 summarizes the number 
of absorption components from the original absorber catalogue that 
meet successive selection criteria. All components in our final sample 
are detected at >3o significance, and none are strongly saturated (i.e. 
the central optical depth of C Iv 41550 is t < 2). 


2.4 Grouping components 


The C iv absorber statistics can be calculated using either individual 
Voigt profile components of the recovered C Iv absorbers, or systems 


4z = 5.46 is the path-length-weighted mean redshift of the survey when 
quasar proximity zones are included. 
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Table 1. The number of absorption components from the orig- 
inal absorber catalogue that meet successive selection criteria, 
as described in Section 2.3. 


Criterion Number of components 
Civ 727 
+ Primary Sample 559 
+ Intervening 507 
+ 13.2 < logN < 15.0 273 


that group all components lying within a chosen velocity interval Av. 
The choice of whether to use components or systems directly impacts 
the measured number density and the slope of the column density 
distribution function (e.g. Boksenberg & Sargent 2015). The latter 
translates into a minor impact on &cıy because it is measured over a 
specific log N range (see Section 3.3). 

In this work, we choose to examine the statistics of individual 
components. The number density and column density distribution 
function of individual components are strongly resolution-dependent 
because components separated by less than the spectral resolution 
may be blended. In Section 3, we directly compare our C IV absorber 
statistics with those of lower redshift absorbers identified in much 
higher resolution spectra (R ~ 50000). To enable consistent com- 
parison of all data sets, we merge components with Av < 50kms^! 
(this threshold is commonly used because it is the coarsest spec- 
tral resolution for typical absorption line measurements). Merged 
components are combined into a single ‘system’ defined by the 
total column density and column-density-weighted mean redshift 
of the constituent components (Songaila 2001; D'Odorico et al. 
2010). We emphasize that the merging is not intended to group 
all physically associated absorption components, but rather provides 
a means to consistently compare data sets that have significantly 
different spectral resolutions. The 273 C Iv components in our final 
sample comprise a total of 260 systems. 


2.5 Comparison with previous catalogues 


Catalogues of CIV absorption have been previously published for 
12/42 of the E-XQR-30 quasars, based on either archival spectra 
or lower quality spectra of the XQR-30 quasars (Codoreanu et al. 
2018; Meyer et al. 2019; D'Odorico et al. 2013, 2022). A detailed 
comparison with these catalogues is given in section 5.4 of Paper I. 
In brief, we recover 107/122 (88 per cent) of the previously reported 
CTV systems. 4/15 (27 per cent) of the missed systems were found to 
be better explained by transitions of other ions at different redshifts. 
Our comprehensive simultaneous search for low and high-ionization 
systems allows us to more robustly identify the origin(s) of individual 
absorption lines compared to searches focused on individual ions. 
One system (7 percent) was rejected by the checkers due to strong 
skyline contamination. The remaining 10/15 (60 percent) missed 
systems have column densities near or below the completeness limit 
of our data and were missed due to low S/N or blending with much 
stronger absorption features. D’Odorico et al. (2013) performed line 
identification visually and D'Odorico et al. (2022) filtered candidates 
manually rather than using an automated algorithm, and this likely 
explains why they were able to detect some weak systems that were 
missed in our catalogue. However, the improved quality of the E- 
XQR-30 spectra enabled the detection of approximately 50 new C Iv 
systems along the lines of sight to the same quasars, increasing the 
number of known systems in these spectra by ~30 per cent. 
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Figure 2. The absorption path-length for intervening C Iv absorbers covered 
by the E-XQR-30 survey over the two redshift bins shown in the left-hand 
panel of Fig. 1 (blue and red pentagons). Black markers show the path-lengths 
covered by other surveys of CIV absorption in the literature. The literature 
samples differ significantly from one another in the column density thresholds 
above which they are complete. 


2.6 Absorption path 


The first step in calculating the incidence of C Iv absorbers over a 
given redshift interval is to compute the total absorption search path 
covered by the survey. The absorption distance to an object at a 
redshift of z is given by: 


2 2 
XO) = za (AmC +z’ + 2a)” (1) 


m 

(Bahcall & Peebles 1969). The absorption search path covered by a 
single sightline for which C Iv is accessible over the range zı € z < 
z2 is given by AX = X(z1) — X(z»). If the redshift search interval is 
broken up by masked or BAL regions, then AX should be computed 
for each subinterval individually. The total absorption path covered 
by our sample of 42 quasar spectra over a redshift range zı € z < z2 
can be calculated by taking the redshift search intervals shown by the 
black horizontal lines in the left-hand panel of Fig. 1, and summing 
AX for all subintervals that lie within that redshift range. The redshift 
search intervals can be reproduced using the data in tables B1, B2, 
and B4 of Paper I which list the adopted quasar redshifts, the masked 
wavelength regions, and the BAL regions (originally published by 
Bischetti et al. 2022 and Bischetti et al. submitted), respectively. The 
absorption path-length values quoted in this paper can be recovered 
using the PYTHON code? and data provided with the metal absorber 
catalogue (see footnote 1). 

Fig. 2 compares the absorption path-length for intervening CIV 
absorbers covered by the E-XQR-30 survey to those of previous 
surveys of CIV absorption in the literature. To aid in comparison, 
we divide our sample into two redshift bins (4.317 < z < 5.235 and 
5.236 < z < 6.338) that have approximately equal path-lengths, and 
align well with the redshift intervals covered by previous surveys. 
Our sample has a total absorption path of AX = 105.9 in both 
redshift bins; a factor of ~2.6 (71.4) larger than the previous largest 
data set at z = 5.2 (4.3 S z S 5.2). This data set represents the largest 


The routine can be used to calculate the absorption path for the public E- 
XQR-30 metal absorber catalogue over a given redshift interval for a range 
of ions/transitions and any value of Qm. 
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Table 2. Completeness-corrected C Iv absorber statistics measured after grouping components with Av < 50km s^. The statistics are reported in two 
column density ranges and four redshift bins chosen to cover approximately equal absorption path-lengths. The Qcyy values given in this table are 
calculated using equation (10). We provide two sets of error estimates for each of dn/dX and Q2c y. The SL96 approximation for ó€2c y is calculated using 
equation (11). In each error column, the two listed values represent the lower and upper Io error. 


ó(dn/dX) ó(dn/dX) Qew Nc (x 1078) BQeq (x 1078) 
z range <z> AX # dn/dX (bootstrap) (poisson) (x 1078) (bootstrap) (SL96) 
13.2 < log (N/cm?) < 15.0 
4.317 <z < 4.915 4.72 52.93 115 2.69 0.39, 0.33 0.25, 0.27 2.63 0.66, 0.59 0.43 
4.915 <z < 5.236 5.07 52.94 74 1.73 0.22, 0.23 0.20, 0.23 1.37 0.27, 0.28 0.19 
5.236 <z < 5.590 5.40 52.94 42 0.96 0.16, 0.17 0.15, 0.17 0.94 0.16, 0.21 0.18 
5.590 <z < 6.339 5.7] 52.96 29 0.79 0.23, 0.13 0.15, 0.18 0.55 0.22, 0.16 0.14 
13.8 < log (N/cm?) < 15.0 
4.317 <z < 4.915 4.72 52.93 24 0.48 0.19, 0.14 0.10, 0.12 1.52 0.63, 0.57 0.42 
4.915 <z < 5.236 5.07 52.94 13 0.26 0.08, 0.09 0.07, 0.09 0.67 0.21, 0.26 0.17 
5.236 <z < 5.590 5.40 52.94 12 0.24 0.06, 0.07 0.07, 0.09 0.56 0.17, 0.18 0.17 
5.590 <z « 6.339 5.7] | 52.96 8 0.16 0.08, 0.06 0.06, 0.08 0.32 0.19, 0.14 0.14 
homogeneous set of CIV absorbers at z = 5 measured at S/N = 10 Age (Gyr) 
with spectral resolution FWHM x30 km s™!. 4 3 15 12 1 
13.2 < logN < 15.0 
3 Civ LINE STATISTICS ai ¢ 
* 
3.1 Number density : a 
The number density dn/dX is defined as the number of absorbers per a T 4 
unit comoving path-length. We investigate how the number density I y 
of Civ absorbers evolves over the probed redshift range by dividing 924 + Fin ME Correction y 
the sample into four redshift bins, listed in Table 2. The bins were € D'Odorico+10 (re-calculated) a 
"M ; * D'Odorico+22 (re-calculated) * 
chosen to divide the survey path-length approximately equally and 1d v Codoreanu-18 (logN > 12.5) : J 
provide measurements of dn/dX with a time resolution of ~100 Myr, A E 13 (logN 7 13.4) B. 
enabling us to robustly constrain the rate of evolution in the number sds Finlator--20 (C boosted by 1.75) — 
density of C iv absorbers over 4.3 S z S 6.3. 0 r T T T — 
PPS 2 3 4 5 6 
The measured number density is strongly dependent on the z 


considered column density range. We calculate dn/dX over the entire 
column density range of our sample (13.2 « log N « 15.0). To enable 
completeness corrections to be applied, the absorbers in each redshift 
bin are divided into column density bins with a width of 0.1 dex (28 
bins in total), and dn/dX is computed as follows: 


28 


dn "T n(log N;) 
ax 2 = » C(z,log Nj) @) 


Here, n(log N;) represents the number of absorbers in the ith log N 
bin and C(z, log Nj) is the sample completeness at the relevant z 
and log N. The latter is computed using the arctan fits provided in 
Paper I (see equation 3 and table 2 therein) which are shown in the 
right-hand panel of Fig. 1. We adopt the lower redshift completeness 
curve for all measurements at z < 5.235 and the higher redshift curve 
for all measurements at z > 5.235. We note that very similar results 
are obtained when calculating the completeness correction factors 
directly from the mock absorber catalogues (following the method 
described in Paper I) over the exact redshift range of each bin. 

We estimate the errors on dn/dX both assuming Poisson statistics 
and using bootstrap resampling. For each redshift bin, the bootstrap 
error is calculated by selecting quasar sightlines for which the absorp- 
tion path-length over that redshift bin is non-zero, randomly sampling 
t times from this set of t sightlines (with replacement), measuring 
dn/dX, repeating this 200 times, and computing the standard deviation 
of these measurements. We also report Poisson errors based on the 
lo upper and lower bounds on 7 calculated using equations (9) and 
(14) from Gehrels (1986), respectively. Our measurements are listed 


Figure 3. Number density of C IV absorbers over the column density range 
13.2 € log N < 15.0 in four redshift bins (orange markers). The open and solid 
markers show measurements before and after correcting for completeness, 
respectively. Horizontal error bars indicate the widths of the redshift bins and 
vertical error bars are calculated using bootstrapping. Black markers show 
measurements from the literature with minimum column densities indicated 
in the legend. The grey curves show predictions from the Finlator et al. 
(2020) cosmological simulation over the same redshift and column density 
ranges, before (solid) and after (dotted) rescaling the carbon abundances (see 
Section 3.4). 


in Table 2. The bootstrap errors are typically larger than the Poisson 
errors because the bootstrap sampling accounts for factors such as 
variation in S/N (and therefore completeness) between spectra and 
cosmic variance. We adopt the boostrap errors in the subsequent 
analysis. 

The number density measurements are shown in Fig. 3. The 
open and filled orange markers show the values before and after 
correcting for completeness, respectively, and the black markers 
show measurements from other studies in the literature. The statistics 
for the D'Odorico et al. (2010) and D'Odorico et al. (2022) samples 
have been re-calculated to match the log N ranges used in this paper. 
The errors on z and dn/dX indicate the widths of the redshift bins and 
the bootstrap errors, respectively. 

We find that the number density of C Iv absorbers decreases by a 
factor of 3.4 + 0.9 over the ~300 Myr interval between z ~ 4.7 and 
~5.8. Interestingly, ~91 per cent of this drop occurs over the 200 Myr 
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Figure 4. Left: The Civ CDDF in two redshift bins centred at z ~ 4.9 (blue) and z ~ 5.6 (red). Within each redshift bin, the absorbers are split into 11 log N 
bins each with a width of 0.3 dex spanning 11.7 < log N « 15.0. The open and solid markers show measurements before and after correcting for completeness, 
respectively. We do not compute completeness-corrected values in the range where the survey is incomplete (log N < 13.2, brown shaded region). Vertical error 
bars represent the uncertainties assuming Poisson statistics. Upper limits are reported for bins where no absorbers are detected. We use a maximum likelihood 
approach to fit a power law distribution to the unbinned absorber data at each redshift. The red and blue-dashed lines and shaded regions show the best fit and the 
associated 68 per cent confidence interval at each redshift, respectively. Grey curves show the CDDFs of simulated absorbers from the Finlator et al. (2020) data 
set over the same two redshift ranges as our measurements. Right: lo (solid), 2e (dashed), and 30 (dotted) confidence intervals of the joint posterior probability 
distribution functions for o (slope) and B (normalization) for both redshift intervals. We measure a strong decrease in the normalization of the CDDF (0.38 dex) 
but no significant evolution in the slope, demonstrating that the decline in C Iv absorption occurs approximately uniformly across 13.2 < log N < 15.0. 


interval between z ~ 4.7 and ~5.4. Our results are qualitatively 
consistent with the findings of D'Odorico et al. (2022) who report 
a sharp drop in the C Iv number density between z ~ 4.7 and ~5.3, 
and Codoreanu et al. (2018) who measured a significant decrease in 
the C IV number density between z ~ 4.8 and ~5.7. 

The large number of C IV absorbers used in this study enables us 
to robustly confirm that the decline in CIV absorption over 4.3 S z 
S 6.3 is much more rapid than what is observed over 1.5 < z < 4. A 
factor of 3.4 decrease in dn/dX over a period of 300 Myr corresponds 
to approximately a factor of 1.5 decrease in dn/dX every 100 Myr. A 
similar decrease in number density (a factor of 1.5 + 0.3) is measured 
over 1.8 S z € 3.4 from the D'Odorico et al. (2010) sample, but this 
covers an order of magnitude more cosmic time (1.75 Gyr) than the 
E-XQR-30 sample. The number density we measure in the z — 4.7 
bin is very similar to the value measured at z ~ 3.4, suggesting that the 
rapid decline in dn/dX likely commences at z ~ 5. Our results clearly 
show that the number density of C Iv absorbers decreases slowly over 
1 S z $5 before declining rapidly towards higher redshifts. 


3.2 Column density distribution function (CDDF) 


The column density distribution function (CDDF) is a fundamental 
property of an absorber population, similar to the galaxy luminosity 
function. It is defined as the number of absorption systems per unit 
column density per unit redshift absorption path dX, and is denoted 
as f(N). A population of absorbers with constant distributions of 
physical cross-sections and comoving space densities illuminated by 
a non-evolving ionizing spectrum will have the same CDDF at all 
redshifts. 

Measuring the CDDF at different redshifts is important because 
it can reveal whether the observed drop in CIV number density 
is dominated by absorbers in a particular column density range. 
If only the weakest (strongest) absorbers are suppressed at higher 
redshifts, then the CDDF would be flatter (steeper) at earlier times. 
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To accurately constrain the slope of the CDDF, we require as many 
absorbers per redshift interval as possible. Therefore, we characterize 
the CDDF for the two redshift intervals shown in Fig. 2 (centred at z 
~ 4.9 and ~5.6). 

To visualize the CDDF at each redshift, we first measure f(N) 
without completeness corrections over the whole column density 
range covered by the sample of primary intervening C Iv absorbers 
(11.7 < log N < 15.0). Within each redshift bin, the absorbers are 
split into 11 log N bins with a width of 0.3 dex. We also measure 
completeness-corrected values for bins lying above the 50 percent 
completeness limit (log N > 13.2) by further splitting each log N bin 
into three sub-bins of width 0.1 dex and computing f(N) as follows: 


1 y n(log N; ;) 
AN;AX C(z, log Ni, j) 


j=1 


(Ni, z) = logio (3) 


Here, AN; is the column density range (in linear space) covered 
by the ith log N bin. When no absorbers are detected in a given bin, 
we report an upper limit based on the Poisson upper bound for n 
(Gehrels 1986). 

The f(N) measurements are shown in the left-hand panel of Fig. 4 
and listed in Table 3. The values before and after correcting for 
completeness are shown as open and solid markers, respectively. 
These binned measurements are not used when fitting the CDDF (as 
explained below), so we report only Poisson errors. The CDDF rises 
in the region where the sample is incomplete (brown shaded region) 
before declining towards larger column densities. As expected, the 
completeness corrections have a relatively minor impact on f(N) for 
bins above the completeness limit. 

The binned measurements help to visualize the slope and normal- 
ization of the CDDF in the two redshift bins. However, to minimize 
biases associated with binning, we fit the CDDF using the binning- 
independent maximum likelihood approach outlined in Bosman et al. 
(2017). We parametrize the CDDF as a power law distribution with 
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Table 3. Completeness-corrected measurements of f(N) in two redshift 
intervals. Within each redshift interval, the absorbers are split into log N 
bins of width 0.3 dex. When no absorbers are detected we report an upper 
limit based on the Poisson upper bound on the number of absorbers. 


log N range # logf(N) dlog AN) (poisson) 
4.317 < z < 5.235 (<z> = 4.914) 
13.20-13.50 92 —13.13 0.06,0.04 
13.50-13.80 60 —13.69 0.06,0.06 
13.80-14.10 19 —14.51 0.11,0.11 
14.10-14.40 12 —15.03 0.15,0.14 
14.40-14.70 4 —15.81 0.28,0.25 
14.70-15.00 2 —16.42 0.45,0.36 
5.236 < z < 6.338 (<z> = 5.589) 
13.20-13.50 3| —13.52 0.12,0.05 
13.50-13.80 20 —14.16 0.11,0.10 
13.80-14.10 14 —14.65 0.13,0.13 
14.10-14.40 4 —15.51 0.28,0.25 
14.40-14.70 2 —16.12 0.45,0.36 
14.70-15.00 0 < —16.46 - 
the form 
N —a 
Nes). " 
FN) (=) (4) 
We adopt Ny = 10/?9* for consistency with previous works 


(D'Odorico et al. 2013; Codoreanu et al. 2018), but note that our 
conclusions remain the same regardless of the adopted value. The 
likehood function is written as follows: 


L(B, a) = P,(n|B, o) | | vito). (5) 
The first term P,(n|B, œ) represents the Poisson probability of 
observing a total of n absorbers across 13.2 « log N < 15.0 given 
a CDDF with parameters (B, o), for which the expected number of 
absorbers is given by 


1015.0 
C(z,N 
Nexpectea( B, a) = BAX x 10% f G N) 


dN. 6 
1013.2 N” (6) 


The second term P;(N;|o) represents the Poisson probability of 
observing a system with column density N; when B is scaled such 
that P,(n|B, œ) = 1. We calculate this probability by comparing the 
number of observed and expected absorbers in bins of width 0.1 dex, 
because the mean uncertainty on the individual log N measurements 
is 0.09 dex. However, we note that adopting a bin size of 0.01 dex 
produces consistent results. 

We perform the fitting using the Affine Invariant Markov Chain 
Monte Carlo (MCMC) Ensemble Sampler EMCEE (Foreman-Mackey 
et al. 2013). We adopt uniform priors on o and B over the ranges 
—20 < logio B < —10 and 0 < o < 4. We run the MCMC with 300 
walkers, 100 burn-in steps, and 500 run steps. The initial parameter 
estimates for the walkers are distributed uniformly over the ranges 
—13 x log Bio < —17 and 0 <a <4. 

The right-hand panel of Fig. 4 shows the lo, 2c, and 30 contours 
of the joint posterior probability distribution function for œ and B, 
and the best-fitting values and lo errors on both parameters are listed 
in Table 4. The dashed lines and shaded regions in the left-hand panel 
of Fig. 4 show the best-fit power law distributions and the 68 per cent 
confidence intervals, respectively. 

We find that the normalization of the CDDF (B) decreases strongly 
over cosmic time, changing by AB = —0.38 + 0.06 dex (a factor 
of 2.4 + 0.3) between z ~ 4.9 and ~5.6. This is consistent 
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with the observed drop in number density over the same redshift 
interval. Our findings are in agreement with previous studies which 
report consistent decreases in normalization over a similar redshift 
range (Codoreanu et al. 2018 measure AB = —0.47 + 0.20dex 
and D'Odorico et al. 2013 report a factor of ~2-3 decrease in 
normalization). However, the large size of the E-XQR-30 sample 
enables us to measure this evolution at a significance of >5o for the 
first time. The Codoreanu et al. (2018) sample includes 41 systems 
over four sightlines and the D'Odorico et al. (2013) sample includes 
102 systems over six sightlines at 4.35 — z < 6.2, compared to the 
current sample of 260 systems over 42 sightlines. 

In contrast, there is no evidence for evolution in the slope of 
the CDDF over 4.3 S z S 6.3. We measure power law slopes of œ 
= 2.197995 atz 4.9 ando. — 2.18 +0.15 at z ~ 5.6, corresponding 
to a change in slope of Aw = —0.01 + 0.17. This is consistent with 
previous works which also found Ao ~ 0 over a similar redshift range 
(e.g. Aa = —0.18 + 0.36; D'Odoricoetal. 2013, Aœ = 0.47 €: 0.47; 
Codoreanu et al. 2018). However, our measurements provide a more 
stringent constraint on Aq due to the significant increase in sample 
size. 

We note that the individual o values measured in this work are 
larger than the values reported by D'Odorico et al. (2013) and 
Codoreanu et al. (2018) which range from 1.44 to 1.96. This likely 
stems from a combination of factors. First, there are some differences 
in adopted methodology between the studies. Codoreanu et al. (2018) 
group components with Av « 500 km s^! which is expected to result 
in a flatter CDDF than the 50 km s^! grouping used in this work (e.g. 
Boksenberg & Sargent 2015). If we group components with Av 
< 500kms~!, we measure o ~ 1.8-2.0, in better agreement with 
the values published by Codoreanu et al. (2018). D'Odorico et al. 
(2013) performed least squared fitting on the binned measurements 
in contrast to the unbinned maximum likelihood approach adopted 
in this work. Finally, the samples used in the previous studies had 
similar completeness limits to our sample but very few C Iv absorbers 
with log N > 13.8 (14 in D'Odorico et al. 2013 and 8 in Codoreanu 
et al. 2018 over both redshift intervals combined). The current data 
set includes 54 such systems, increasing the dynamic range in column 
density over which reliable measurements can be made and thereby 
improving the constraint on o. 

Put together, our results provide the first robust indication that the 
observed decrease in the C IV dn/dX is not dominated by absorbers 
in a particular column density range but occurs approximately 
uniformly for absorbers spanning 13.2 < log N « 15.0. 


3.3 Cosmic mass density (Qc) 


The Civ cosmic mass density (Gc,,) is defined as the CIV mass 
per unit comoving Mpc, pcıy, divided by the critical density of the 
Universe at z ~ 0, pai. The comoving mass density of CIV can be 
obtained by integrating the CDDF: 


Homc 


[e 


Pcl = 


J vronan. (7) 


where Ho is the Hubble constant and mc is the atomic mass of carbon. 
The C IV cosmic mass density is then given by: 


Qe = pctv/ Perit (8) 


Storrie-Lombardi, McMahon & Irwin (1996, hereafter SL96) pro- 
posed that the integral of the CDDF can be approximated as UN/AX, 
yielding the following expression: 
Homc XN 

C Perit AX Í 


Qew X 


(9) 
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Table 4. Parameters of the power law distributions (equation 4) fit to the unbinned C Iv absorber statistics using the maximum 
likelihood approach described in Section 3.2, and Qcjy values computed from the integral of the CDDF (equation 7). The errors on 
a and f represent the 1o confidence intervals measured from the posterior probability distribution functions shown in the right-hand 
panel of Fig. 4. The errors on Qcıy represent the 68 per cent confidence interval of values obtained by calculating Qc, for every 


realization of o and £. 


Qciy(x 1075) 
z range <z> log B a 13.2 < log N/cm? < 15.0 13.8 < log N/cm? < 15.0 
4.317 < z < 5.236 4.914 —13.76+0:03 2.19008 1.99*02 15 a 
5.236 < z < 6.339 5.589 —14.14+0:04 24802 0,815058 0.471045 


Many previous studies of C Iv absorption at z ~ 5-6 have used this 
approximation to measure €2c;, due to small sample sizes (e.g. Ryan- 
Weber et al. 2009; D'Odorico et al. 2013; Díaz et al. 2016; Meyer 
et al. 2019; D'Odorico et al. 2022). We present measurements of 
Qew obtained using both the integral of the CDDF and the SL96 
approximation. 

The measured Qc,y is strongly dependent on the adopted log N 
range. Many previous studies report CIV absorber statistics over 
13.8 < log N < 15.0, but our sample is 750 per cent complete all the 
way down to log N= 13.2. Therefore, we present Qc;y measurements 
summed over both 13.2 < log N < 15.0 and 13.8 < log N < 15.0. 

The Qe, values calculated from the integral of the CDDF are 
already completeness-corrected because the completeness is taken 
into account when fitting the power law distribution (see equation 6). 
The errors on these measurements are computed by calculating Qc y 
for each (B, o) pair in the MCMC joint posterior distribution and 
taking the standard deviation of these values. The measurements 
obtained from equation (7) and the corresponding MCMC errors are 
listed in Table 4. 

We use the SL96 approximation to compute Qc, for the four 
redshift bins listed in Table 2. Once again, to enable completeness 
corrections, the absorbers in each redshift bin are divided into column 
density bins with a width of 0.1 dex and Qc yy is calculated as follows: 


MNT Homc > Mootal,i 
CN Cpa AX — C(z, log Ni)’ 


i 


(10) 


Here, Mota; is the total column density of all absorbers in the ith 
log N bin. The uncertainty on Qc: is typically computed using the 
following approximation for the fractional variance, also proposed 


by SL96: 
_ X(N?) 


(Sen) 

Qeiv (ENy' 
This approximation is widely adopted in the literature but is known to 
underestimate the true errors by up to a factor of 1.5 (e.g. D'Odorico 
et al. 2010). As previously discussed, bootstrapping provides a 
more accurate estimate of the errors for sufficiently large samples 
because it accounts for factors such as variations in completeness 
between spectra. Table 2 lists the Qc; measurements calculated 
using equation (10) along with errors calculated using both the 
SL96 approximation (equation 11) and bootstrapping. We adopt the 
bootstrap errors in our analysis. 

The results are shown in Fig. 5. The left-hand and right-hand 
panels show Qc ;y measurements summed over 13.2 < logN < 
15.0 and 13.8 < logN < 15.0, respectively. The measurements in 
the left-hand panel are based on a larger number of absorbers (see 
Table 2) and therefore have smaller associated errors than those in 
the right-hand panel. The orange markers with and without black 
borders show the measurements computed from the integral of the 
CDDF (equation 7) and the SL96 approximation (equation 10), 


(11) 
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respectively. The two methods produce consistent results. In each 
panel, the black markers show literature measurements for which the 
minimum log N is consistent with our adopted value to within 0.2 dex. 
The literature measurements have been scaled to match our adopted 
cosmology. 

Our results are consistent with previously published values, but the 
larger sample size enables us to measure Qc: in a larger number of 
redshift bins and therefore track the time evolution more finely than 
previous studies. We find that Qc, declines steeply but smoothly 
between z ~ 4.7 and ~5.8, decreasing by a factor of 4.8 + 2.0 
over this interval of only 300 Myr for absorbers with 13.2 < log N < 
15.0. (The values measured over 13.8 < log N < 15.0 yield consistent 
results but at a lower statistical significance due to the smaller sample 
size). 

Combining our results with measurements at 1.5 « z « 4 from 
D'Odorico et al. (2010), we find that Qc, evolves much more rapidly 
at z = 4.3 than at lower redshifts. Qc decreases by a factor of 
3.5 + 1.3 over the 1.75 Gyr interval from z ~ 1.8 to ~3.4; comparable 
to the decline between z ~ 4.7 and ~5.8 which occurs over a much 
shorter time interval. The Qc ,y values measured at z ~ 4.7 and ~3.4 
are consistent with one another, suggesting that, similar to dn/dX, the 
rapid decline in Qcy commences at z = 5. 

Our results are consistent with many previous works which also 
suggested that Qc, remains approximately constant over 1.5 Sz <5 
before declining rapidly towards higher redshifts (e.g. Songaila 1997, 
2001, 2005; Becker et al. 2009; Ryan-Weber et al. 2009; D'Odorico 
et al. 2010; Simcoe 2011; D'Odorico et al. 2013; Boksenberg & 
Sargent 2015; Díaz et al. 2016; Bosman et al. 2017; Codoreanu 
et al. 2018; Meyer et al. 2019; D'Odorico et al. 2022). However, 
the improved measurements presented in this paper constrain the 
rate and timing of the accelerated decline in Qc, to much higher 
statistical significance than was previously possible. 


3.4 Comparison to predictions from the Technicolor Dawn 
simulation 


We compare our results to predictions from the Technicolor Dawn 
cosmological simulation (described in detail in Finlator et al. 2020). 
The simulation is run using a custom version of GADGET-3 and 
covers a volume of 15 h~! Mpc at a mass resolution of 2.6 x 10? Mo. 
Metal enrichment is modelled using the supernova yields of Nomoto 
et al. (2006) scaled to the assumed Kroupa (2001) IMF with a 
hypernova fraction of 50 percent. Stellar feedback is implemented 
using momentum Kicks tuned to match the Muratov et al. (2015) 
outflow mass-loading factor scaling, but with a slight reduction 
in normalization to better match the observed galaxy stellar mass 
and UV luminosity functions. The outflow velocities are taken from 
Muratov et al. (2015) with two boost factors applied to produce larger 
outflow velocities that are more consistent with observed values 
(see also Davé, Thompson & Hopkins 2016). The quasar radiation 


€z0z 1snBny 9z uo sen Áq $626904/v€/L/LzS/ermie/seuujfuoo: dno'oruepeoe//:sdyu Woy pepeojuwoq 


1.2 1 0.8 


PX 
n - E 13.2 < logN < 15.0 


2 € EXQR-30 
QG 108] 6 EXQR3O (from CDDF) 
* D'Odorico+22 (re-calculated) 
€ D'Odorico+10 (re-calculated) 
« Simcoe-11 (logN > 13.4) 
x  D'Odorico-13 (logN > 13.4) 
m Bosman+17 (logN > 13.4) 
--- Delayed Enrichment (t = 610 Myr) 
—— Finlator+20 
senes Finlator+20 (C boosted by 1.75) 


10-9 T r r r 


2 3 4 5 6 7 


Zz 


E-XQR-30: C iv absorption over 4.3 <z < 6.3 323 


Age (Gyr) 
1.5 1.2 1 


1077 1 1 — L 1 1 
13.8 < logN < 15.0 


2 0-8 
O 10-734 
Gq 
» Ryan-Weber+09 (logN > 13.8) 
* Meyer+19 (logN > 13.8) 
+ Diaz+16 (logN > 13.8) 
x D’Odorico+13 (logN > 13.8) 
10-9 r : r 


2 3 4 


Figure 5. Redshift evolution of Qcıy measured over 13.2 < log N < 15.0 (left) and 13.8 < log N < 15.0 (right). The orange markers with and without black 
borders represent measurements made using the integral of the CDDF (equation 7) and the SL96 approximation (equation 10), respectively. Orange-dashed 
curves show the best-fit delayed enrichment model for each set of measurements (for the Madau & Dickinson 2014 fit to the cosmic star formation rate density, 
as described in Section 4.1). Other plotting symbols and errors are the same as in Fig. 3. Literature measurements of €2c y have been scaled to match our adopted 
cosmology. The large E-XQR-30 data set enables us to robustly measure C2c y at intervals of ~ 100 Myr, revealing a smooth but steep decline from z ~ 4.7 to 


5.8. 


field is assumed to be spatially uniform while the inhomogeneous 
galaxy UV background is modelled by solving the multifrequency 
radiative transfer equation for galaxies on a Cartesian grid. A mock 
absorber catalogue was generated by passing sightlines through the 
simulation volume. The metallicity of the gas particles is traced self- 
consisently in the simulation, while the ionization state is computed 
using ionization equilibrium calculations (see Finlator et al. 2015). 

The mock absorber catalogue has previously been compared 
with measured absorber statistics at 3 S z S 5 (Hasan et al. 2020) 
and 5 S z S 6.3 (Finlator et al. 2020). At 3 S z S 5, Technicolor 
Dawn roughly reproduces the observed statistics of C Iv absorbers 
with W > 0.3 À but significantly over-produces weaker absorbers, 
perhaps indicating that the simulated outflows transport metals too 
far from galaxies (Hasan et al. 2020). At z = 5, the simulation 
is able to reproduce the observed CDDF of Sirv absorbers but 
produces a factor of ~3 too few CIV absorbers at all column 
densities (Finlator et al. 2020). This tension is also seen in other 
cosmological simulations (e.g. Keating et al. 2016; Rahmati et al. 
2016). Finlator et al. (2020) argue that this discrepancy cannot be 
solved by increasing the strength of the UV background, the star 
formation efficiency or the overall CGM metallicity because all 
of these would lead to discrepancies in other galaxy or absorber 
properties. One potential avenue to reduce the disrepancy could be 
to relax the assumption that H II regions are ionization bounded. 
If a non-negligible fraction of HII regions are density bounded, 
then Lyman continuum photons would be able to escape through 
attenuated channels, hardening the emergent ionizing spectrum. 

An alternative solution could be to artificially increase the yield 
of carbon relative to silicon. The carbon yields are impacted by 
uncertainties on the supernova yields and the initial mass function, 
which may be more top-heavy in early low-metallicity galaxies (e.g. 
Nomoto et al. 2006; Kulkarni et al. 2013). Finlator et al. (2020) 
measured mean C Il/Si I rest-frame equivalent width ratios of 0.54 
(0.77) for simulated z ~ 6 absorbers with Si 1 41260 (C n 1.1334) 
rest-frame equivalent widths exceeding 0.05 A. These values are 22— 
75 per cent lower than the mean ratio of 0.94 measured for observed 
systems at 5.75 < z < 6.25 in the Becker et al. (2019) sample. The 


mean C 1/Si 1 equivalent width ratio for absorbers at 5.75 < z < 6.4in 
the E-XQR-30 catalogue (Paper I) is 1.08, supporting the hypothesis 
that the relative carbon abundances of the simulated absorbers in 
the Finlator et al. (2020) sample could be underestimated. Finlator 
et al. (2020) found that increasing the carbon abundances of the 
simulated absorbers by a factor of 1.75 brings their number densities 
into agreement with the observed statistics. 

We compare the C Iv absorber statistics predicted by the Finlator 
et al. (2020) simulations to the observed statistics in Figs 3, 4, 
and 5. The solid and dotted grey curves represent the predictions 
before and after rescaling the carbon abundances by a factor of 1.75, 
respectively. The curves were calculated using simulated absorbers in 
the same column density intervals and redshift ranges as the observed 
systems. Prior to rescaling the carbon abundances, the simulation 
underproduces C Iv absorbers at all column densities. After rescaling 
the carbon abundances, the predicted number densities are in good 
agreement with our measurements (Fig. 3). The rescaled CDDFs 
are marginally consistent with the measured statistics near the 
completeness limit of our sample (log N — 13.2), but stronger 
systems are increasingly underproduced moving towards higher 
column densities (Fig. 4). Hasan et al. (2020) similarly find that the 
simulation underproduces strong C Iv absorbers at 3 S z < 5, perhaps 
due to the limited simulation volume or the absence of quasar-driven 
spatial fluctuations in the hard UV background. As a direct result, the 
rescaled Qcyy predictions fall a factor of 1.2-2.0 (0.9-1.90) below 
the observed €2c,, computed over 13.2 < log N < 15.0, and a factor 
of 1.8-4.8 (1.1-2.30) below the observed Qc; computed over 13.8 
x log N < 15.0 (Fig. 5). We note that boosting the carbon abundances 
leads to the overprediction of Qc, at z < 4, whereas the original 
predictions are in agreement with the measured values. 

Overall, the discrepancy between the measured and predicted 
normalization of the Civ CDDF supports the hypothesis that the 
Finlator et al. (2020) simulation underestimates the carbon yield 
of early stellar populations. However, our robust measurements of 
strong absorber statistics enabled by the large E-XQR-30 sample size 
additionaly reveal that the simulation assembles strong systems too 
slowly. 
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4 DISCUSSION 


The statistical power of the combined sample of 42 quasars con- 
sidered in this work has enabled us to robustly characterize the 
evolution of dn/dX, Qc,,, and the slope and normalization of the 
Civ CDDF over 4.3 S z S 6.3. Our results conclusively show that 
the incidence of C IV absorbers decreases rapidly over this redshift 
range. We additionally show that the slope of the C 1v CDDF does 
not evolve significantly (Aw = —0.01 + 0.17), suggesting that the 
suppression of CIV absorption occurs relatively uniformly across 
13.2 x logN « 15.0. In the following sections, we use these new 
measurements to investigate what fraction of the observed Qeyy 
evolution could plausibly be driven by steady enrichment of the 
CGM/IGM by outflows (Section 4.1), and by changes in the average 
ionization state of the absorbing gas (Section 4.2). 


4.1 Chemical enrichment by outflows 


4.1.1 Toy model 


We first explore the possibility that the rapid evolution in cy 
could trace early enrichment by galactic winds. Outflows transport 
metals liberated by supernovae and stellar winds from galaxies to the 
CGM/IGM, steadily increasing the metal content of galaxy haloes 
over cosmic time. Our sample probes C Iv absorption at a time when 
the Universe is $1.5 Gyr old; too young for a significant fraction of 
ejected metals to have been re-accreted onto galaxies (which occurs 
over a time-scale of —1—3 Gyr; Oppenheimer et al. 2009; Brook 
et al. 2014; Übler et al. 2014). Furthermore, the majority of galaxies 
ejecting metals at these redshifts have very shallow potential wells, 
increasing the likelihood that the metals will escape completely (e.g. 
Ferrara 2008; Pallottini et al. 2014). Therefore, it is reasonable to 
assume that the majority of metals that have been ejected remain 
outside of galaxies. 

Following Madau & Dickinson (2014, hereafter MD14), we 
assume that the comoving mass density of metals produced in 
and ejected from stars up to a given redshift z is approximately 
proportional integral of the cosmic star formation rate density (z^): 


IH(z) , dz 
PZ, ejected(Z) = yd = nl Wz ) dr (12) 
0 


Here 1 — R is the fraction of the initial stellar mass that is retained 
in stars, y is the average mass of heavy elements created and ejected 
into the CGM/IGM per unit mass in those stars, and fz) is the age 
of the Universe at redshift z. The mass density of ejected carbon can 
be computed by multiplying equation (12) by the mass fraction of 
metals in carbon, Ac: 


Ac Pz, ejected (Z) (13) 


Finally, the expected mass density of C Iv in the ejected material can 
be obtained by multiplying equation (13) by the fraction of carbon 
which is triply ionized: 


pc, ejected(Z) = 


Civ 
aPC, ejected(Z)- (14) 

Crotal 

In this section, we aim to quantify the expected evolution in Qeyy 

originating from changes in the carbon content alone, so we assume 

that C IV/Ciou is independent of redshift. 

In order to compute the predicted absolute evolution of Qc 
over cosmic time, it is necessary to substitute in values for y, R, 
Ac, and CIV/Ciota. However, the stellar population properties and 
CTV ionization fraction are poorly constrained at these redshifts. We 
circumvent the uncertainties in these parameters by investigating 


pcw(z) ~ 
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whether the predicted rate of change in the mean metal content is 
sufficient to explain the observed rate of evolution of Qcjy, assuming 
that y, R, and Ac do not vary significantly over the ~ 1.25 Gyr between 
the beginning of the Universe at z ~ 4.7. Specifically, we compare 
the predicted and observed fractional growth of cy between two 
redshifts zı and z2: 


Adi D 


Dci G1) " VG) 
pew) -— Tuc. ?|£ |ar- 


We emphasize that the calculation of xci, does not require us to 
adopt a specific value for the carbon yield or the fraction of produced 
carbon that is ejected into the CGM/IGM. However, it does require 
the extrapolation of V (z^) to z > 8, where observational constraints 
are currently relatively limited. In the future, large surveys with 
JWST will enable much more robust measurements of y/(z') at these 
redshifts. We account for the uncertainty on v (z^) by performing our 
calculations using two parametrizations that differ significantly at 
z = 8. The MD14 parametrization has previously been used to 
calculate the cosmic mean metal content as a function of redshift 
(e.g. Péroux & Howk 2020). However, constraints from the growing 
number of very high redshift galaxies observed with HST and JWST 
suggest that it overpredicts y at z = 8 (e.g. Mason, Trenti & Treu 
2015; Oesch et al. 2018; Donnan et al. 2023; Harikane et al. 2022b, 
a). Early JWST results additionally suggest that distant galaxies are 
relatively blue, making it unlikely that a significant fraction of star- 
formation is obscured (e.g. Castellano et al. 2022; Cullen et al. 
2023; Ferrara, Pallottini & Dayal 2022; Finkelstein et al. 2022; 
Nanayakkara et al. 2022; Santini et al. 2022). We therefore also 
consider a second y(z’) parametrization published by Harikane et al. 
(2022b, hereafter H22b) (their equation 60) which assumes a constant 
star-formation efficiency at z > 10. This parametrization captures the 
observed rapid decline at z > 6 but the scaling is too steep and under- 
estimates the measured y at z = 10 (e.g. Bouwens et al. 2022a, b; 
Harikane et al. 2022a). As a result, the xc; values implied by the 
MD14 and H22b parametrizations may bracket the true fractional 
growth of Qc, as a result of chemical enrichment. 

Based on the MD14 (H22b) parametrization for w(z’), the ex- 
pected fractional growth in the mean carbon mass density between 
z = 5.77 and 4.72 (given by xcw (5.77, 4.72)) is 2.1 (2.8). The H22b 
parametrization predicts a steeper rise in y(z’) at high redshift than 
the MD14 model, and therefore also predicts a faster growth in the 
cosmic mean metal content. However in both cases, the expected 
change in pc is significantly smaller than the observed factor of 
4.8 evolution in Qeyy. 

This discrepancy may be resolved by accounting for the time delay 
between the formation of metals in stars and the deposition of these 
metals into the CGM/IGM where they are observed as C Iv absorbers. 
Low-mass (S8 Mo) stars eject significant quantities of metals during 
the asymptotic giant branch (AGB) phase, with a typical time delay 
of hundreds of Myr (e.g. Oppenheimer & Davé 2008; Kramer, 
Haiman & Madau 2011). It may then take a similar amount of time 
for outflows to transport these metals from the interstellar medium 
(ISM) to the sites of C Iv absorption which are typically located on 
scales of tens to hundreds of kpc (e.g. García et al. 2017a; Díaz et al. 
2021; Kashino et al. 2022). To account for this, we introduce a delay 
time T into equation (15) such that the metal mass density of the C 
IV-absorbing gas at a redshift z is assumed to be proportional to the 


dt 


F3 


Xew(Z1, 22) = (15) 


integral of w(z’) evaluated at tg; — T: 
foe" we) | | ae f 
Xew(Z1, Z2, T) = : (16) 
R 2)— v n) & | dr 
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Figure 6. Left: Measured xcıy (see equation 15) for absorbers with 13.2 < log N < 15.0 (orange pentagons), and the best-fit delayed enrichment curves for 
200 bootstrap realizations of the measurements for the MD14 parametrization of the cosmic star formation rate density y(z’) (coloured curves). The colour of 
each curve indicates the delay time of the corresponding best-fitting model. The delayed enrichment model assumes that Q2c ;y reflects the metal content of the 
CGM/IGM and scales as the integral of y(z’) with some delay time t (see equation 16). Right: Distribution of best-fit delay time values for the 200 realizations, 
for both the MD14 (green) and H22b (blue) parametrizations of y(z’). The dashed lines indicate the median values of the distributions (250 and 605 Myr). The 
MD14 parametrization has been shown to overpredict the star formation rate density at z = 8 whilst the H22b predicts too little star formation, suggesting that 


^ 


the delay time required to explain the Qcıy evolution by chemical enrichment alone likely lies between these two values. 


4.1.2 Fitting and physical interpretation 


We investigate which delay time value best explains the observed 
Qe evolution for each w(z’) parametrization as follows. We 
consider 90 different delay times spanning t = 0-900 Myr with 
At = 10 Myr. For each delay time t, we compute xc y (5.77, z) for 
the four values of z corresponding to the centres of our four redshift 
bins (given in Table 2), and then calculate the total X^ difference 
between the four measured xc, values and the predicted values. The 
delay time yielding the smallest x? value is then selected as the best- 
fiting value. The uncertainty on the best-fit delay time is estimated 
by repeating this calculation for all 200 bootstrap estimates of Qc 
(see Section 3.3). 

The xci, measurements are shown as orange pentagons in the left- 
hand panel of Fig. 6. To illustrate the range of delay time values 
allowed by the bootstrap measurements, we overplot the best-fit 
delayed enrichment curves for the MD14 parametrization, where the 
colour represents the delay time of the best-fit model according to the 
colourbar to the right of the plot. In the right-hand panel, we show the 
distribution of best-fit delay times for the 200 bootstrap iterations, for 
both the MD14 (green) and H22b (blue) models. The upper x-axis 
shows the implied formation redshifts of the metals observed at z 
~ 5.8. These plots are based on the Qcıy values measured over log N 
= 13.2-15.0, but consistent results are obtained for Qc,;y measured 
overlog N — 13.8-15.0 (with larger errors due to the smaller sample 
size). 

We find that the observed Qc ,, evolution can be fully explained by 
gradual enrichment with a median best-fit delay time of 605 Myr for 
the MD14 parametrization, or 250 Myr for the H22b parametrization. 
As discussed in Section 4.1.1, these two parametrizations appear to 
bracket the true cosmic SFR density at z = 10, and therefore it seems 
likely that the true delay time required to explain the €2c;, evolution 
by chemical enrichment alone lies somewhere between these two 
values; i.e. 250 Myr S t S 605 Myr. 

In this scenario, the metals observed in C Iv absorbers at z ^ 5.8 
would have been formed at redshifts of z ~ 7-12; between ~ 100 and 
7-600 Myr after the expected first light of galaxies at z ~ 15-20 (e.g. 


Hashimoto et al. 2018; Labbe et al. 2022; Naidu et al. 2022; Yan et al. 
2022). The abundance patterns of the highest redshift metal absorber 
at z = 6.84 suggest that early populations of Population I-like stars 
were already present at z = 10 (e.g. Simcoe et al. 2020). However, 
it is unclear whether it is physically reasonable to invoke such long 
delay times. 

We consider two sources of delay: the time taken for stars to 
return their carbon to the ISM, r,4, and the time taken for outflows 
to transport metals from the ISM to the CGM/IGM where they 
are observed as CIV, Tou. Kramer et al. (2011) investigated the 
distribution of delay times between the formation of stars and the 
return of carbon to the ISM for a single age stellar population, 
including contributions from both Type II supernovae and AGB stars. 
They found that the mean delay time for a Chabrier (2003) IMF is 
Tret = 160 Myr.° Assuming that the outflow delay time Tow is given 
by tou & t — Tret, Our results imply that Tou ~ 90—445 Myr. 

The lower end of this Touw range is consistent with the predictions 
of Oppenheimer et al. (2009) who find that simulated C IV absorbers 
at z ~ 6 are typically observed 30-300 Myr after being launched in a 
wind. Touw can also be expressed as the quotient of the characteristic 
distance from galaxies at which CIV absorbers are detected (the 
impact parameter D) and the outflow velocity Vout; 1.6. Tout = D/Vout. 
This relationship can be used to infer the characteristic impact 
parameter of CIv absorbers for a given delay time. Observations 
of C Il emission and rest-frame UV absorption suggest that star- 
forming galaxies at z ~ 5-6 have typical outflow velocities of 300— 
500 km s^! (e.g. Gallerani et al. 2018; Sugahara et al. 2019; Ginolfi 
et al. 2020; Pizzati et al. 2020). Over travel times of 90-445 Myr, such 
outflows would transport metals to typical distances of 30-230 kpc. 
These values can be considered as upper limits because outflows may 
decelerate as they propagate (e.g. Nelson et al. 2019). Oppenheimer 
et al. (2009) find that the typical velocities of winds associated with 
z ~ 6, CIV absorbers, averaged over their lifetimes, are ^50 per cent 


$250 Myr for a Kroupa (2001) IMF 
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of their launch velocities. Accounting for this deceleration suggests 
that the ejected metals would more likely reach impact parameters 
of 15-115 kpc. 

The lower half of this impact parameter range is consistent with 
predictions from simulations and with observed values. Oppenheimer 
et al. (2009) predict that most C IV absorbers at z ^- 6 trace metals on 
their first journey into the IGM, making it simpler to link absorbers 
with galaxies at this time compared to later epochs. They find that 
most simulated C Iv absorbers at this redshift lie ~5-50 kpc from 
their host galaxies. Observationally linking galaxies to absorbers 
at these redshifts is very challenging because dwarf galaxies are 
expected to play a significant role in the enrichment of the CGM/IGM 
(e.g. García et al. 2017a; Diaz et al. 2021). However, several studies 
have identified galaxies within 750 kpc of strong C Iv absorbers. Cai 
et al. (2017) find a Lyo emitter (LAE) 42 kpc from a C Iv absorber 
with log N > 14.0 at z ^ 5.74. Diaz et al. (2021) find that 2/8 of the 
faint LAEs associated with C Iv absorbers with log N > 13.5 at 4.9 <z 
< 5.7 have impact parameters of <50 kpc. Codoreanu et al. (2018) 
find that C Iv systems with log N > 14 are preferentially associated 
with low-ionization counterparts and may therefore trace multiphase 
gas at the interface between the CGM and the IGM. Some strong C Iv 
absorbers do not appear to have any galaxies in this impact parameter 
range but instead show one or more galaxies at larger distances, 
suggesting that the absorbing gas could be part of an extended 
structure connecting galaxies (e.g. Díaz et al. 2021). However, the 
presence of intrinsically faint or dust-obscured galaxies below the 
detection limit at smaller impact parameters cannot be ruled out 
(e.g. Neeleman et al. 2020). Deeper searches for galaxies associated 
with C IV absorbers are required to obtain better constraints on their 
typical impact parameters. 

Together, these findings suggest that the observed rapid evolution 
in Qcw over 4.3 S z S 6.3 could plausibly be explained by 
gradual chemical enrichment with a delay time of 250-605 Myr, 
comprising ~160 Myr for the release of carbon from AGB stars and 
7-90-445 Myr for the transport of carbon in outflows. Assuming 
lifetime-averaged outflow velocities of —150—250 km s^!, the ex- 
pected typical impact parameter for C Iv absorbers would be —15— 
115 kpc, which is consistent with observed and predicted values. 

The orange-dashed lines in Fig. 5 show the predicted Qc 
evolution across 4 S z S 7 for the best-fit MD14 model. (By design, 
the curve representing the best-fit H22b model is almost the same). 
If the carbon enrichment of the CGM progresses in this manner, then 
CIV absorption should become increasingly rare at z > 6, which is 
consistent with the scarcity of detections at these redshifts. Cooper 
et al. (2019) and Simcoe et al. (2020) find no evidence of C Iv in the 
earliest known intervening metal absorber at z — 6.84, and Bosman 
et al. (2017) measured an upper limit on Qc at z ~ 6.5 from a single 
z ~ 7 quasar sightline (see left-hand panel of Fig. 5). 

However, we note that this model cannot explain the observed evo- 
lution in Qeyy atz < 4. The model predicts that cy should continue 
to increase rapidly at z S 5; in tension with the measurements which 
show that Qc,, changes very little between z ~ 4.7 and ~3.4. As 
mentioned in Section 4.1.1, the toy model assumes that the current 
metal content of the CGM/IGM represents the accumulation of all 
metals ejected from galaxies in the past, and this assumption starts 
to break down when the Universe becomes old enough that a non- 
negligible fraction of the ejected metals have been re-accreted onto 
galaxies. However, it seems unlikely that metal recycling would lead 
to such a pronounced change in the growth rate of Qc y between z ~ 5 
and ~4. In Section 4.2.3, we discuss the possibility that this transition 
may be more naturally explained by changes in the ionization state 
of the absorbing gas. 
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4.1.3 Caveats 


The calculations presented here are based on a toy model with 
many associated uncertainties, some of which have already been 
discussed. The w(z) measurements on which the MD14 and H22b 
parametrizations are based only include contributions from galaxies 
with Muy < —17, corresponding to SFR > 0.3 Mo yr^!. However, 
fainter galaxies have shallower gravitational potentials and may 
therefore contribute substantially to the enrichment of the CGM/IGM 
(e.g. Garcia et al. 2017a; Diaz et al. 2021). 

We assume that the carbon yields of stars remain constant between 
the beginning of the Universe and z ~ 4.7. However, studies of metal- 
poor stars in the Milky Way have revealed a significant population of 
carbon enhanced metal poor (CEMP) stars (see Beers & Christlieb 
2005 for a review), and there is evidence for enhanced [C/O] ratios in 
the most metal-poor stars and metal absorbers (e.g. Cooke, Pettini & 
Steidel 2017; Bafiados et al. 2019). If stellar populations at z ~ 6 have 
enhanced carbon yields compared to stellar populations at z ~ 5, the 
carbon content would grow more slowly than the integral of y(z’), 
requiring larger delay times to explain the observed Qc,y evolution 
which could become unrealistic given the age of the Universe at these 
redshifts (^41 Gyr at z ~ 6). 


4.2 Ionization state of the CGM 


We also consider the possibility that the rapid evolution in Qc, could 
be driven by a change in the ionization state of the CGM/IGM. It is 
unclear whether gas in the haloes of z ~ 5-6 galaxies is primarily 
ionized by the UV background or by radiation from the host galaxy. 
As discussed earlier, simulations often underpredict the incidence of 
CTV absorbers at z Z 5, which could be an indication that galaxies 
self-ionize their haloes (e.g. Keating et al. 2016). If the halo gas 
is primarily ionized by the local radiation field, then the ionization 
state of the CGM/IGM would be expected to correlate with the typical 
strength and hardness of this radiation field. 

However, observational evidence suggests that the UV background 
isthe primary source of ionization in galaxy haloes. Lau, Prochaska & 
Hennawi (2016) found that the ionization parameter of the gas 
around quasar host galaxies at z ~ 2-3 increases towards larger 
radii, suggesting that the halo gas is primarily ionized by the UV 
background. Meyer et al. (2019) found evidence for excess Lyo 
transmission in the vicinity of C Iv absorbers at z ~ 5—6, suggesting 
that these absorbers trace local enhancements in the UV background 
(see also Finlator et al. 2016). If the UV background is the dominant 
source of ionization, then the properties of metal absorbers may be 
significantly impacted by the reionization of hydrogen and/or He IL. 


4.2.1 Heu reionization 


The energy required to ionize Hell (54.4 eV) lies between the 
energies to ionize CIN and C v (47.9 and 64.5 eV, respectively), 
and therefore He II reionization may impact the ionization states of 
carbon absorbers. He H reionization is expected to end at around z 
~ 3 (e.g. Davies & Furlanetto 2014; Upton Sanderbeck, D' Aloisio & 
McQuinn 2016; Worseck et al. 2019; Makan et al. 2021). Songaila 
(1998) reported a jump in the Si1V/C Iv ratio at z ~ 3 that they 
attributed to a change in the UVB associated with the end of Hell 
reionization, but subsequent studies did not reproduce this finding 
(e.g. Kim, Cristiani & D'Odorico 2002; Boksenberg, Sargent & 
Rauch 2003; Cooksey et al. 2011; Boksenberg & Sargent 2015). 
The onset of He II reionization is uncertain but may have occurred 
around z ^ 4.5, when the temperature of the IGM has been observed 
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to rise (e.g Becker et al. 2011; Boera et al. 2014; Walther et al. 2019; 
Gaikwad et al. 2021). This is intriguingly similar to the redshift where 
we observe a transition in the rate of evolution of Qcyy. It could be 
that the properties of CIV absorbers are impacted by changes in 
the relative contribution of QSOs to the UV background and the 
mean free path of Hell ionizing photons during He II reionization. 
However, the incidence rate of Sirv (which is not sensitive to the 
He ll break in the UV background) follows the same trends as C Iv 
(D’ Odorico et al. 2022). Furthermore, the C Iv incidence rate changes 
very little over z ~ 3-4.5 when the Hell filling factor of the IGM 
must be changing rapidly. For these reasons, we suggest that He II 
reionization is likely not a key factor driving the observed evolution 
in Qc. 


4.2.2 Hydrogen reionization 


The end of hydrogen reionization is expected to be marked by rapid 
changes in the amplitude of the UV background (e.g. Kulkarni et al. 
2019; Keating et al. 2020; Garaldi et al. 2022; Lewis et al. 2022), 
driving a transition of metals from primarily neutral or low-ionization 
environments to more highly ionized environments once reionization 
is complete (e.g. Becker et al. 2019). There are multiple independent 
lines of evidence suggesting that signatures of reionization persist 
below z ~ 6 (e.g. Zhu et al. 2021; Bosman et al. 2022), indicating 
that this effect could be significant over the redshift range probed 
in this study. Indeed, there is growing observational evidence for a 
shift in the typical ionization environments of metal absorbers at z 
~ 5-6. This work confirms previous findings that the incidence of 
Civ absorbers declines rapidly at z = 4.9, and similar behaviour 
is observed for Siiv (e.g. D'Odorico et al. 2022). Examining the 
evolution of low-ionization ions can help to deduce whether the 
decline in high-ionization absorbers is primarily driven by ionization 
or metal content. The incidence of weak Mg II absorbers remains 
approximately constant over the probed redshift range (e.g. Bosman 
et al. 2017; Chen et al. 2017; Codoreanu et al. 2017; Zou et al. 2021), 
and the incidence of weak O1 absorbers increases over 4.1 <z < 
6.5 (e.g. Becker et al. 2019). The cosmic mean metallicity decreases 
towards higher redshift, so this upturn can only be explained by an 
increase in the fraction of oxygen observed as O I (e.g. Doughty & 
Finlator 2019) and/or an increase in the relative oxygen abundance 
(see e.g. Welsh et al. 2022 and references therein). 

It is interesting to consider whether both the upturn in weak (W 
< 0.2 À rest-frame) O I absorbers at z = 5.7, and the decline in 
CTV absorption at z = 5 can be explained by a single population of 
absorbers transitioning from higher ionization environments at z ~ 5 
to lower ionization environments at z ~ 6. In this scenario, if weak 
C1 absorbers follow the same trend as weak O I absorbers, then 
we would expect the decline in C Iv absorption to be dominated by 
systems with W < 0.2 A, which corresponds to an apparent column 
density of log Na (C Iv) S 13.7 (assuming the absorption lies on the 
linear part of the curve of growth; Spitzer 1978). However, we have 
shown that the slope of the C rv CDDF does not vary as a function 
of redshift (see Section 3.2 and Fig. 4), indicating that the strong 
and weak C Iv absorbers decrease in incidence at similar rates. We 
caution that the measured number densities of strong (W > 0.2 À) Or 
absorbers are relatively uncertain due to small number statistics and 
do not preclude the possibility of an upturn at z > 5.7 (e.g. Becker 
et al. 2019). The interpretation is additionally complicated by the 
fact that O I probes dense gas whereas C IV probes diffuse gas which 
in lower ionization environments may manifest as C III A977 Å (e. g. 
Cooper et al. 2019) which falls in the saturated Lyo forest. 
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Figure 7. Bar graph illustrating the contribution of CIV and CH to the 
cosmic mass density of carbon in each redshift bin. Vertical lines illustrate 
the measurement errors and are centred on the path-length-weighted mean 
redshift of each bin with small offsets added for clarity. We are unable to 
measure Qc y in the two lowest redshift bins due to the saturation of the Lyo 
forest. The evolving balance between €2c; and Qe across the two highest 
redshift bins suggests that changes in the UV background driven by hydrogen 
reionization may contribute to the decline in Qc y. 


4.2.3 Constraints from measurements of Qc 


Changes in the ionization state of the absorbing gas can be examined 
more directly by measuring the column density ratios of different ions 
ofthe same species (e.g. log M(C 11)/log N(C Iv)) which are insensitive 
to changes in metallicity and/or abundance patterns over cosmic time. 
Previous studies have found evidence for an increase in the C 1/C Iv 
ratio at z = 5.7 (e.g. D'Odorico et al. 2013; Cooper et al. 2019; 
Simcoe et al. 2020), suggesting that some portion of the observed 
Qcıy evolution could be the result of carbon transitioning to lower 
ionization states. 

We use the E-XQR-30 metal absorber catalogue to compute Qcn 
in the same redshift bins as our Qc, measurements and investigate 
whether the decrease in Qc, is accompanied by an increase in the 
frequency of carbon in lower ionization states. We select intervening 
CH absorbers from the primary absorber sample using the method 
described in Section 2.3 and consider systems in the same log N 
range as the CIV absorbers.’ We do not present measurements of 
S2c for the two lowest redshift bins where C II is inaccessible due 
to the saturation of the Lyo forest. In a forthcoming paper, we will 
investigate the use of Mg II as a proxy for C II (e.g. Cooper et al. 2019) 
to place similar constraints on the ionization state of the absorbing 
gas in the lower redshift bins (Sebastian et al., in prep). 

Fig. 7 illustrates the contributions of C Iv and C II (when available) 
to the cosmic mass density of carbon over the four redshift bins used 
in this paper. Over the two highest redshift bins, the decrease in Qcıv 
is completely balanced by an increase in Qcy. At z ~ 5.8, we measure 
Qc = (8.417222) x 107? and Qcnpciy = (1.397023) x 1078, and 
at z ~ 5.4, we measure cu = (4.1 1252) x 107° and Qeniciv 
= (1:35*025) x 1075. This suggests that changes in the typical 
ionization state of carbon as a result of reionization may play a 
significant role in driving the evolution of Qc at the highest redshifts 
probed by our sample. 


7We note that the sample completeness exceeds 50 per cent for C 1 absorbers 
with 13.2 < log N « 15; see Paper I. 
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On the other hand, changes in Qc, cannot balance the evolution 
in Qcyy at lower redshifts. The combined Qcyicry measured at z 
~ 5.4 is smaller than the €2c;y measured at z ^ 4.7, indicating 
that Qcy+crv must increase over this redshift range. We note that 
9] percent (81 percent) of the observed evolution in dn/dX (Qc) 
occurs over this interval. These results suggest that a change in the 
typical ionization state of metal absorbers is likely not the only factor 
driving the evolution in Qcyy over 4.3 S z S 6.3. 

The main limitation of this analysis is that we are only able to 
Observe two ionization states of carbon. It is unclear what fraction 
of the missing CIV might be present in C III which is often the 
dominant ionization state at z — 1 (e.g. Lehner et al. 2018), but 
cannot be observed for any of the objects in our study because 
Cir A977 À falls in the saturated Lya forest. Many simulations 
predict that the overall fraction of carbon in C Iv changes relatively 
little over the probed redshift interval (e.g. Oppenheimer et al. 2009; 
Cen & Chisari 2011; Finlator et al. 2015; Rahmati et al. 2016; García 
et al. 2017b). However, some of these simulations were run using a 
uniform UV background which has been shown to be inconsistent 
with observations of the Lyo forest at z > 5.3 (Bosman et al. 2022). 
The simulations that do model an inhomogenous reionization are 
computationally limited to volumes smaller than the spatial scales 
over which the Lyo forest opacities have been observed to fluctuate 
(Becker et al. 2015b) due to the need to simultaneously model 
galaxy evolution processes on orders-of-magnitude smaller scales. It 
therefore remains unclear whether or not cosmological simulations 
predict a strong evolution in the fraction of carbon observed as C Iv 
across this redshift range. 


4.2.4 Physical intepretation 


Changes in the ionization state of the absorbing gas may offer a 
natural explanation for the observed two-stage evolution of Qew, 
characterized by a rapid rise in Qc, over 6 = z 2 5 followed by 
a much more gradual increase over 5 = z 2 1.5. Bosman et al. 
(2022) found that spatial fluctuations in the UV background and/or 
the properties of the IGM are required to explain the scatter in the Lyo 
optical depth down to z = 5.3. If these fluctuations are associated with 
the end of hydrogen reionization, then we would expect to observe 
a large change in the amplitude of the UV background at similar 
redshifts (e.g. Kulkarni et al. 2019; Keating et al. 2020; Garaldi 
et al. 2022; Lewis et al. 2022), which in turn could contribute to 
the observed rapid evolution of Qc, in the early Universe. The 
flatter evolution in Qcıy observed at z S 5 would then reflect the 
more gradual changes in the UV background following the end of 
reionization. 

To illustrate how a reduction in the amplitude of the UV back- 
ground could lead to a decline in Qcyw, we performed simple 
photoionization modelling with CLOUDY (Ferland et al. 2017) to 
estimate the gas densities probed by different carbon ions at z ~ 5 and 
~6. The gas temperature is fixed at T = 10^ K which is well below 
the collisional ionization temperature for C IV. We adopt the 2011 
update of the Faucher-Giguére et al. (2009) ionizing spectrum, with 
a density-dependent correction for gas self-shielding as described 
in Keating et al. (2016). The amplitude of the ionizing spectrum 
is re-scaled to match the HI photoionization rate measurements 
from Calverley et al. (2011).* The shape and amplitude of the UV 


8PH = 107125 s7! atz = 5 and Tyi = 1071284 s7! at z = 6. This evolution is 
quite steep but is reproduced in some simulations of a late-ending reionization 
(e.g. Lewis et al. 2022). 
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Figure 8. Illustration of the possible ranges of gas densities probed by 
different ionization states of carbon at z — 5 and 6. The curves are the outputs 
of CLOUDY photoionization models. We fix the gas temperature to 10* K and 
adopt the 2011 update of the Faucher-Giguére et al. (2009) ionizing spectrum 
with a density-dependent correction for gas self-shielding as described in 
Keating et al. (2016). The ionizing spectrum is re-scaled at each redshift to 
match the H I photoionization rates measured by Calverley et al. (2011). 


background are difficult to measure at these redshifts and may also 
vary spatially (e.g. Davies & Furlanetto 2016). We emphasize that 
the adopted values are used only to illustrate how changes in the 
strength of the UV background translate to changes in the distribution 
of carbon between ionization states. 

The results are shown in Fig. 8. As the amplitude of the UV 
background decreases (moving from z ~ 5 to ~6), CIV absorption 
(shown in red) becomes increasingly sensitive to diffuse gas which 
may be located further from the galaxies responsible for driving 
outflows and may therefore be less metal-enriched, helping to explain 
the observed decline in €2c y. The highest density material probed by 
CIV absorption at z ~ 5 is predominately traced by C mM (green) at 
z ~ 6, highlighting the importance of considering the unobservable 
Cm absorption when investigating the redshift evolution of Qc. 
It is also clear that the range of densities probed by C II absorption 
(orange) extends towards lower values at z ^ 6 than at ~5. These 
effects may partially explain the observed decline in the Qey/Qcen 
ratio over the two highest redshift bins of our sample (Fig. 7). We note 
that the shape of the ionizing spectrum may evolve simultaneously 
with the amplitude between z ~ 6 and ~5 (e.g. Finlator et al. 2018). 
This effect is neglected here, but would likely enhance the trends 
seen in Fig. 8. 

It is not straightforward to determine how changes in the UV 
background would impact the slope of the Civ CDDF. Keating 
et al. (2016) found that the relative number of strong and weak 
CTV absorbers is sensitive to both the shape and amplitude of the UV 
background. However, they did not compute their simulated CDDF 
and quantify the changes in its slope. Finlator et al. (2016) found 
that any changes in the CDDF slope for different UV background 
models would be too small to detect given current measurement 
uncertainties, even for an extreme model where the UV background 
is completely dominated by QSOs. More theoretical investigation is 
therefore required to understand the impact of the evolving amplitude 
and shape of the UV background in the context of a late-ending 
reionization, and whether this would produce a detectable signature 
in the slope of the CDDF. 
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5 SUMMARY 


We use a sample of 260 C Iv absorbers at 4.3 S z S 6.3 to examine 
the rate and physical origin of the rapid decline in the cosmic mass 
density of CIV (Qc) at z = 5. The absorbers were drawn from 
the catalogue published by Davies et al. (2023) who performed 
a systematic search for metal absorbers in high S/N XSHOOTER 
spectra of 42 z ~ 6 quasars, of which 30 are drawn from the ESO 
large program XQR-30. The sample contains a factor of ~3 more C 
IV absorbers at z > 5 than any previous compendium in the literature, 
enabling us to make the most robust measurements to date of C Iv 
absorber statistics in the high-redshift Universe. 

We measure the number density (dn/dX) and cosmic mass density 
(Qcıw) of CIV absorption in four redshift bins. The sample is 
>50 percent complete for log N > 13.2, and our analysis is focused 
on absorbers in the column density range 13.2 < logN < 15.0. 
Both dn/dX and Qc, decline rapidly but smoothly over the probed 
redshift range, with dn/dX decreasing by a factor of 3.4 + 0.9 and 
Cc, decreasing by a factor of 4.8 + 2.0 over this ~300 Myr period. 
This decline is much more rapid than what is observed over 1.5 < z 
« 4.0 for C Iv absorbers from the D'Odorico et al. (2010) sample. 
Interestingly, the statistics of C IV absorbers at z ~ 4.7 are comparable 
to those at z ~ 3.4. Our results are consistent with previous findings 
that the incidence of strong C Iv absorbers is approximately constant 
over 1.5 S z S 5 and declines rapidly towards higher redshifts. 
However, the larger sample size used in this work enables us to 
measure dn/dX and Qc, in smaller redshift intervals, and therefore 
constrain the rate and timing of the decline to higher precision than 
previous studies. 

We investigate whether the observed evolution is dominated by 
absorbers in a particular column density range by fitting the CIV 
CDDF in two redshift bins centred at z ~ 5.6 and ~4.9. The 
normalization decreases by a factor of 2.4 + 0.3, consistent with 
the observed number density evolution, but there is no evidence for a 
change in slope between the two redshift bins (Ao = —0.01 + 0.17). 
Our work provides the first robust indication that the decline in C IV 
absorption occurs approximately uniformly across 13.2 < log N < 
15.0 and is not preferentially driven by strong or weak absorbers. 

We use our improved measurements to examine whether the rapid 
evolution in Qc, could reflect a change in the carbon content of 
the CGM/IGM driven by the fast production of metals in the high- 
redshift Universe. Assuming that the mean metal content of the 
CGM/IGM scales with the integral of the cosmic star formation rate 
density, it would be expected to increase by a factor of <2.8 between 
z ~ 5.8 and —4.7, which is less than 60 per cent of the observed Qc y 
evolution. If we additionally account for the time taken to deposit 
carbon into the ISM and eject it into galaxy haloes (7250-605 Myr), 
we find that chemical evolution could plausibly explain all of the 
observed evolution in Qc, over 4.3 S z S 6.3. However, this model 
cannot easily explain why 2c, declines much more slowly over 1.5 
SzS5thanatz> 5. 

Finally, we explore whether the cy evolution could be driven 
by a change in the fraction of carbon observed as C IV by comparing 
the relative evolution of Qc and Qcy. Examining two ions of the 
same element removes uncertainties associated with variations in 
metallicity and/or abundance patterns over cosmic time. Over the two 
highest redshift bins (z ~ 5.4 and ~5.8), we find that the decrease in 
Qcıv is completely balanced by an increase in Qc, suggesting that 
the decline in Qc y over this redshift range may trace a decrease in the 
typical ionization state of carbon as the IGM becomes more neutral. 
In contrast, the combined £2c,, + Qc measured at z ~ 5.4 is smaller 
than Qceyy alone at z ~ 4.7. This could indicate that a change in the 
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typical ionization state of metal absorbers is not the dominant factor 
driving the evolution in Qc, over 4.3 S z S 6.3. However, we are 
unable to measure the contribution from C m which is inaccessible 
at these redshifts. We note that rapid changes in the amplitude and/or 
hardness of the UV background as a result of hydrogen reionization 
may offer a natural explanation for the increased rate of decline in 
Scy at z Z 5. 

Much of this interpretation is based on toy model estimates for 
the gas-phase metallicity and measurements of the Qen/Qcew ratio 
as a proxy for the ionization state of the absorbing gas. In the future, 
photoionization modelling of column density measurements for ions 
covering a wide range of ionization states could be performed to 
obtain probability distributions for the metallicity and ionization 
parameter of the absorbing gas at z ~ 5 and ~6 given different 
assumptions for the strength and shape of the ionizing photon 
background (e.g. Glidden et al. 2016; Cooper et al. 2019; Simcoe 
et al. 2020). Furthermore, carbon emission lines (including CO, [C 1] 
and C irr]) are a complementary probe of carbon in early galaxies. 
Combining emission line measurements with our absorption line 
results would provide further constraints on the evolution of the total 
carbon abundance in the early Universe. 

Overall, the new measurements presented in this paper suggest that 
gradual chemical enrichment by outflows could plausibly explain the 
observed rapid evolution in Qc;y over 4.3 S z S 6.3, without any need 
for changes in the ionization state of the absorbing gas. However, the 
increased incidence of C I1 and O1 absorbers in the early Universe 
(see also e.g. Becker et al. 2019; Cooper et al. 2019) combined with 
the steepening of the decline in Qc, at z 2 5 suggest that reionization 
may also have a significant impact on the evolution of C Iv absorbers 
over cosmic time. 
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